Dynamic spin response of a strongly interacting Fermi gas 



o 

t)JO 

< 



I 



I 

o 
o 



> 

in 

en 

o 

X 



S. Hoinka\ M. Lingham^ M. Delehaye^'^, and C. J. Vale^ 

^ Centre for Atom Optics and Ultrafast Spectroscopy, 

Swinburne University of Technology, Melbourne 3122, Australia 

^ Departement de Physique, Ecole Normale Superieure, 24 rue Lhomond, 75005 Paris, France 

(Dated: February 28, 2013) 

We present an experimental investigation of the dynamic spin response of a strongly interacting 
Fermi gas using Bragg spectroscopy. By varying the detuning of the Bragg lasers, we show that it is 
possible to measure the response in the spin and density channels separately. At low Bragg energies, 
the spin response is suppressed due to pairing, whereas the density response is enhanced. These 
experiments provide the first independent measurements of the spin-parallel and spin-antiparallel 
dynamic and static structure factors which provide insight into the different features of density and 
spin response functions. At high momentum the spin-antiparallel dynamic structure factor displays 
a universal high frequency tail, proportional to uj~^'^ , where huj is the probe energy. 
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Two-component Fermi gases near Feshbach resonances 
provide a well controlled setting to explore many-body 
phenomena in highly correlated quantum systems [T, 
[2]. When the interparticle interactions are sufficiently 
strong, ultracold atomic gases display universal features, 
where macroscopic parameters become independent of 
the microscopic details of the interatomic potential ^- 
[S]. Most studies to date have focussed on static aspects 
of universality |5t410|: however, certain dynamical prop- 
erties can also become universal. Key among these are 
dynamic susceptibilities which describe the way a system 
responds to a perturbation. Recent theoretical work has 
shown that the dynamic structure factor, which is con- 
nected to the imaginary part of susceptibility through the 
fluctuation-dissipation theorem, shows a universal high 
frequency tail both at high momentum |11H13| and low 
momentum where it depends on the frequency dependent 
shear viscosity |14[ [TS] . 



Bragg spectroscopy is a well established tool to mea- 
sure dynamic and static density response functions [16- 
I19| . Previous work on Fermi gases has shown that the 
static structure factor follows a universal law which arises 
from Tan's relation for the density-density correlator [20J. 
Several theoretical studies have also investigated the dy- 
namic spin response [12J [T3] [5T]-|25] and a recent study of 
universal spin transport examined the static spin suscep- 
tibility |24[ [25] , yet the dynamic spin susceptibility has 
not been studied experimentally. 

In this letter we present the first measurements of the 
dynamic spin response of a strongly interacting Fermi 
gas. Two-photon Bragg scattering is used to probe ei- 
ther the spin or density response by appropriate choice 
of the Bragg laser detuning. This allows full characterisa- 
tion of the spin-parallel and spin-antiparallel components 
of the dynamic and static structure factors through the 
application of the /-sum rule [2H1|27]. The spin response 
is suppressed at low energies due to pairing and displays 
a universal high frequency tail, decaying as w~^'^, where 



huj is the probe energy (Bragg frequency) |13j . 

The key to accessing the spin response in two-photon 
scattering experiments is to use Bragg lasers with a dif- 
ferent coupling to each of the two spin states in the mix- 
ture. This can be achieved using spin-fiip Bragg spec- 
troscopy |22l E5] , polarisation sensitive coupling [51] or 
by detuning the Bragg lasers close to resonance. In our 
experiments with ^Li, the first two methods prove chal- 
lenging because of the atomic state configuration at high 
magnetic fields, and therefore we use the third method. 

To understand our measurements we first review the 
atomic level structure of ^Li atoms. Atoms are cooled to 
degeneracy in an equal mixture of the two ground states 
\F = 1/2, mF = ±1/2) labelled |t) and \i). At mag- 
netic fields near the broad Feshbach resonance at 834 
G the electronic and nuclear spins are almost fully de- 
coupled and atoms undergo transitions which preserve 
the nuclear spin projection mi. Both states |t) and \D 
have electronic angular momentum 771,7 = ^1/2 which 
can couple to excited states with mj' = 1/2,-1/2 and 
—3/2. At these magnetic fields the splitting, tu^^, be- 
tween states It) and |4,) is approximately 80 MHz and 
the splitting between excited electronic states is of order 
1.5 GHz. Choosing the polarisation of the Bragg beams 
so that only a~ transitions are possible {mj = —1/2 to 
mji = —3/2), combined with the large Zeeman split- 
ting between excited states, it is straightforward to iso- 
late one, effectively closed, atomic transition from each 
ground state to play any role in the Bragg process. 

In the experiments which follow we work in a regime 
where the Bragg lasers are weak enough to not signif- 
icantly deplete the cloud. We also use a long Bragg 
pulse such that its Fourier width is narrow compared 
to the spectral features being measured. In this limit 
the momentum transferred by the Bragg lasers is propor- 
tional to the imaginary part of the dynamic susceptibility, 
X"(k,a;) =7r[5'(k,a;)-S'(-k,-a;)] HUIII], where ^(k, w) 
is the dynamic structure factor. When the magnitude of 



the Bragg wavevector |k| is large compared to the Fermi 
wavevector kp and the inverse of the de Broghe wave- 
length A^g only the positive term, ^(k, w), contributes 
to the measured response. Thus we can directly obtain 
S{\i,u}) from a Bragg spectrum. In three-dimensional 
Fermi gases 5'(k, w) depends only |k| — k. 

The perturbation introduced by the Bragg lasers can 
be expressed as the sum of two terms, one which couples 
to the total density p{k) ~ p^{k) + Pi{k) and another 
which couples to the total spin S{k). The strength of 
the coupling to each spin state is determined by the light 
shift potential ¥„ — \deg\^ E^ / {hA„) , where A^ is the 
detuning from state a = ^,1., d^g is the electric dipole 
operator for the mj = —1/2 to mji — —3/2 transition 
and E is the electric field amplitude. Choosing the z- 
axis parallel to the Feshbach magnetic field the effective 
Bragg Hamiltonian is given by 
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where /off = 1^ + T^ is an effective total intensity given 
by the sum of the light shifts for each spin state, Beff — 
1^ — y^ is an effective magnetic field given by the dif- 
ferential hght shift |2I], and Sz{k) — p\{k) — p^{k) is 
the z-projection of the total spin. The response in the 
density channel is proportional to the total coupling of 
the Bragg lasers and the spin response is proportional to 
the difference between the couplings to the different spin 
states. As oj^^ ~ 80 MHz, detunings close to resonance 
must be employed to access the spin response. This can 
lead to spontaneous emission and heating for relatively 
modest Bragg laser intensities; however, it is still possible 
to access the linear response regime. 

In a spin-balanced two-component Fermi gas the den- 
sity (spin) dynamic structure factor is defined as ^D] 



SD+.{s-){k,uj) = 2 [S^^{k,uj) ± %(fc,a;)] 
and the components of S{k,uj) are given by 
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where Pa{k, t) is the time-dependent density operator for 
state a \mm\. 

For the case of ^Li at large magnetic field, both the 
It) and 14,) states have essentially identical dipole matrix 
elements d^g for the respective Bragg transitions, and 
therefore the Bragg laser detunings A^ and A^ deter- 
mine the relative strength of the density and spin per- 
turbations. In a spin-balanced Fermi gas S^^{k^ui) = 
S^^{k,u]) and S^\^{k,u]) = S^-\{k,uj) [3^, and recalling 
that S{—k,—uj) -^ for k ^ ^Ft^Jb ^'^ ^^^ frequency 
range we measure, the golden rule [26 gives the momen- 
tum transferred to the atoms by the Bragg perturbation 



to be proportional to 
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From Eq. (4) one can see that using a large detuning 
from resonance, such that A-|- « A^(^ (jj^i) couples pri- 
marily to the density channel as the coefficients in front 
of S'ff and S^^ will be approximately equal. Conversely, 
if At = — A^ then the response will be entirely in the spin 
channel. Choosing detunings between these limits gives a 
variable contribution of both spin and density response. 
The quantities on the right hand side of Eq. (4) each 
satisfy sum rules given by [23i 



N 

UjS-\^{k,Uj)dljJ = -^OJr 

ujS^l{k,Lu)du! = 0, 



(5) 
(6) 



where N is the total number of particles, ujr = hk^ /2m 
is the atomic recoil frequency and m is the atomic mass. 
Thus, the first energy weighted moment of any measured 
response must also satisfy the /-sum rule which allows 
the accurate normalisation of our measured Bragg spec- 
tra 1^ . 

In our experiments we prepare ultracold Fermi gases 
in a balanced mixture of the |t) and jj,) states with 
Na = 2 X 10^ atoms per state as described elsewhere [33] . 
These clouds are loaded into a single beam optical dipole 
trap with trapping frequencies of 97 and 24.5 Hz in the ra- 
dial and axial directions, respectively. Bragg scattering is 
achieved by illuminating the cloud with two beams inter- 
secting at an angle of 84° giving LOr/i^n) = 132 kHz and 
a relative Bragg wavevector of k/kp = 4.5. We use Bragg 
laser detunings of A^ (A^) = 584 (662) MHz to measure 
the density response and A^ (A^) — — (-|-)39 MHz for 
the spin. These detunings, particularly in the spin case, 
are relatively small and can lead to heating via sponta- 
neous emission. For significant intensities this could alter 
the density and spin correlations that determine the re- 
sponse. To establish when this occurs we measure the 
Bragg response for a range of Bragg frequencies in both 
the spin and density channels as a function of the Bragg 
laser intensity as described in the supplemental material. 
Knowing the regimes of linear response, it is possible to 
measure the linear density and spin response functions 
using optimised Bragg laser intensities. 

Spectra are obtained by applying a short Bragg pulse 
to a trapped atom cloud and measuring the momentum 
imparted as a function of the Bragg frequency. Follow- 
ing the Bragg pulse the trap is immediately turned off 
and the atoms are allowed to expand for 500 /is before 
an absorption image of atoms in state |t) is taken. The 
imaging laser frequency is then rapidly switched and a 
second image of atoms in state jt) is taken 850 p& af- 
ter the first. Taking separate images of each spin state 
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Figure 1. Dynamic spin and density structure factors of a two- 
component Fermi gas. The spin (filled squares and circles) 
and density response (dashed and dot dashed lines) where 
measured at l/(kFa) = 0.0 and 1.0, respectively, and the 
dotted blue line is the expected response of a noninteracting 
Fermi gas at T = 0. Solid lines are smoothed curves through 
the spin response data to guide the eye and the responses have 
been normalised using the /-sum rule. 



at different times allows us to measure the differential 
center of mass cloud displacement which is insensitive to 
fluctuations in the trap position. As we use a balanced 
mixture, and with the large coUisional coupling between 
the two spin states, the order in which we image states 
It) and 14,) makes no difference to the measured spectra. 

Figure 1 shows Bragg spectra for both the spin (filled 
squares and circles) and density (dashed and dash-dotted 
lines) at unitarity, l/ikpo) = 0.0, where a is the s-wave 
scattering length, and on the BEC side, \/{kFa) = 1.0. 
Also shown for comparison is the Bragg response of an 
ideal Fermi gas at T = (blue dotted line). These spectra 
were obtained at the lowest temperatures achievable with 
our experiment (^ 0.06 Tp at unitarity). Each spectrum 
is normalised such that the first energy weighted moment, 
J LuS(k,uj)duj, is set to be equal to unity. According to 
the /-sum rule [5^1, this is equivalent to dividing each 
spectrum by the total number of particles multiplied by 
the recoil energy. Therefore by normalising in this way, 
and expressing uj in units of uJr, we obtain the dynamic 
structure factor in units of 1/N |20| . 

Both the spin and density responses are significantly 
different to the ideal gas case indicating the strong cor- 
relations present both in the BEC and unitarity regimes. 
The density spectra SD{k,uj) show a strong response at 
u!r/2 due to the scattering of spin-up/spin-down particle 
pairs, as well as, a broad response at higher frequencies 
in the region of uJr due to the scattering of single atoms 
[TO] . At high momentum the collective mode evolves into 
the pairing feature at uJr/2 |30j . which is clear in both 
spectra, but more prominent on the BEC side due to the 
increased likelihood of finding spin-up/spin-down parti- 
cles at small separation. 

Spin Bragg spectroscopy, on the other hand, is not sen- 



sitive to the collective (paired) mode as the positive and 
negative perturbations on the spin-up and spin-down par- 
ticle, respectively, cancel each other, leaving no nett per- 
turbation on a pair. Therefore the spin response Ss{k, uj) 
is suppressed at ujr/2 and instead only shows free atom 
excitations at higher energies. Ss{k,uj) on the BEC side 
has a lower peak and is biased towards higher frequen- 
cies showing a suppression of the spin susceptibility due 
to the increased energy required to remove atoms from 
bound pairs. 




Figure 2. Spin-parallel and spin-antiparallel components of 
the dynamic structure factor of a strongly interacting Fermi 
gas measured at (a) l/{kFa) = 0.0, and, (b) at l/{kFa) ~ 1.0. 
Red upright triangles are the spin-parallel structure factor Sff 
and blue inverted triangles are the spin anti-parallel response 
Sfi- Solid lines are a guide to the eye. Insets show zoomed 
in plots of the high frequency region where Sfi{k,uj) shows 
a universal tail proportional to uj^^'^. The dotted blue lines 
are power law fits to S-fi{k,ijj > 2ijJr) and the dotted red line 
is the negative of these fits, indicating that Sff{k,ijj) only 
approaches the asymptotic behaviour at higher uu. 

The normalised spectra in Fig. 1 can be combined with 
the known Bragg laser detunings according to Eq. (4) to 
give the spin-parallel {{S^^{k,uj), red triangles} and spin- 
antiparallel {{S^i{k,uj), blue circles} components of the 
dynamic structure factor, plotted in Fig. 2(a) and (b) 
at l/{kFa) = 0.0 and 1.0, respectively. These spectra 
reveal the response of the particle density in one state 
to a perturbation of the particle density in either the 
same or opposite spin states and show quite a complex 
structure. We note that in the absence of interactions, 
the spin-parallel response would be identical to the ideal 
gas response shown in Fig. 1 and the spin-antiparallel 
response would be zero for all uj. 

The features seen in S^-f{k, uj) and S^i{k, ui) can now be 
examined and attributed to different excitations. At fre- 



quencies ranging from to slightly above w^/S, both the 
spin-parallel and spin-antiparallel components are essen- 
tially identical for both interaction strengths. This strik- 
ing result shows that in this frequency range the mea- 
sured response is entirely dominated by the scattering of 
spin-up/spin-down pairs. Each pair contains a spin-up 
and spin-down particle, and, since the response is inde- 
pendent of the combination of spin states being probed, 
the density-density correlations for either combination of 
spin states are the same. 

At higher frequencies, however, the spin-parallel and 
spin-antiparallel response functions become very different 
with S^^{k,uj) showing a significant peak near ujr while 
S^l{k, w) becomes negative for w > Wf The positive peak 
in S^-f{k,u) sa Wr) is primarily due to the autocorrelation 
term, that arises from single particle scattering. This fea- 
ture is not present in S^i{k,u}) that always involves cor- 
relations between different particles. S^i{k,u!) becomes 
negative at high frequencies as required to satisfy the 
sum rule, Eq. (6) which reveals an interesting feature 
of the relative dynamics of spin-up/spin-down particles. 
A high frequency fluctuation with momentum k of the 
spin-down density leads to an out-of-phase fluctuation of 
the spin-up particle density. These anti-correlations sup- 
press the high frequency response in the density channel 
but enhance the high frequency spin response. 

For frequencies much greater than ujr S^^{k,uj) 
{S^l{k, uj)} is predicted to display a universal tail propor- 
tional to -|-{— jl/w^/^, where I is Tan's contact param- 
eter [13]. When measuring the density response alone, 
these two dependencies cancel leaving a universal w"^/^ 
high frequency tail in 5*1) ( A;, w) |11| reminiscent of neu- 
tron scattering measurements on superfluid ''^He [32] . The 
insets of Fig. 2(a) and (b) show zoomed in plots of 
S^^{k,uj) and S-fi{k,uj) at high frequency. Despite the 
small signal we find that the tail of S^i{k,uj) is well 
described by an u}"^/^ dependence. A free power-law 
fit to the data (dotted blue line) for uj > 2uJr yields 
S'^^(fc,a;) = (-0.39±0.11)a; 
%(fc,a;) = (-0.14 ± 0.10) a;-2-7±o 

tent with a;~^/^. 

In the limit that k,uj —}■ cxd, the amplitude of the 
tails should provide a measure of the contact; however, 
our fit coefficients are approximately a factor of three 
larger than expected |lll [IB] . This can be partially at- 
tributed to the fact that the measured response is given 
by the convolution of the dynamic structure factor with 
a squared sine function [TT], [19] which will lift the data; 
however, we estimate that this effect alone would not 
explain the discrepancy. It is more likely that our ex- 
periments are not in the strict fc,w — > cx) limit to obtain 
the contact even though the power law describes the data 
well. We note that S^^{k,uj) only approaches the asymp- 
totic behaviour at higher Bragg frequencies (dotted red 
line) so we do not perform a separate fit to S^^{k, w), nor 
do we resolve the aj~^'^ tail in SnikjUj) in our measured 



2.5±o.3 Qjj i-ijg BEC sijje and 
at unitarity, consis- 



frequency range [11]. 

Finally, the spectra in Fig. 2 can be integrated to pro- 
vide the static structure factors S{k). As these spectra 
are expressed in units of 1/A^, integration over uj gives 
5(fc) directly [20]. We find %(fc) {%(&)} to be 1.02 ± 
0.04{0.17±0.04} at unitarity and 1.03±0.04{0.38±0.04} 
at l/{kFa) = 1.0. At the momentum used here S^^{k) 
is expected to be unity |30| consistent with our findings. 
This also indicates that both the spin and density re- 
sponses provide consistent measures for S^i{k) within our 
uncertainties. The uncertainties in the combined spectra 
are dominated by the uncertainty in the spin response, 
which was obtained at low intensities to minimise spon- 
taneous emission. The situation could be improved by 
using a spin mixture with a larger energy separation such 
as |F = 1/2, mp = 1/2) and \F = 3/2, mp = -3/2) |M] . 

In summary we have shown that Bragg spectroscopy 
with a differential coupling to each spin state can mea- 
sure the dynamic spin and density responses of a strongly 
interacting two-copmonent Fermi gas. We have used this 
technique to extract S^^{k,uj) and S^i{k,uj) separately 
which allows us to associate the different features of the 
density and spin response with the scattering of pairs 
and free atoms. The spin-antiparallel dynamic structure 
factor was seen to display a universal high frequency tail 
proportional to w^^". This work opens the way to mea- 
surements of S^^{k,uj) and 5^^(fc,w) at any momentum, 
which will be particularly important for k < kp where 
5^^(fc) < 1. 
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